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Abstract 

We give sharp estimates for the transition density of the isotropic 
stable Levy process killed when leaving a right circular cone. 



1 Introduction 

Explicit sharp estimates for the Green function of the Laplacian in C > domains 
were given in 1986 by Zhao [53] (see also [351 OS])- Sharp estimates of the Green 
function of Lipschitz domains were given in 2000 by Bogdan 11 . Explicit 
qualitatively sharp estimates for the classical heat kernel in C 1,1 domains were 
established in 2002 by Zhang [62] (see also (63] [32], and [3Ql [54] for further 
extensions). Qualitatively sharp heat kernel estimates in Lipschitz domain were 
given in 2003 by Varopulous [59] . 

The development of the boundary potential theory of the fractional Lapla- 
cian follows an analogous path. Green function estimates were obtained in 1997 
and 1998 by Kulczycki and Chen and Song for C 1 ' 1 domains [46l |29] (see [26l 
Corollary 1.8] for the case of dimension one, see also [33]), and in 2002 by 
Jakubowski for Lipschitz domains [45] (see also [SHUT]). In 2008 Chen, Kim 
and Song gave a sharp and explicit estimate for the heat kernel of the fractional 
Laplacian on C 1,1 domains [26], see ([9]) below. In this note we give an extension 
of the estimate to the right circular cones. We also conjecture, in agreement 
with the results [59] , a likely form of the estimate for a more general class of 
domains: 

pf (x,y) « P x {t d > t)P y (T D > t) Pt {x,y) . (1) 

Here pt (x, y) is the heat kernel of the fractional Laplacian on the whole space R d , 
and P x (te> > t) = J Rd pf(x, y)dy is the survival probability of the corresponding 
isotropic a-stable Levy process. The main result of this paper, Theorem [I] 
asserts that ([1} holds indeed for the right circular cones for all t > 0, x, y € 
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W 1 (see also (|23l) for a more explicit statement). Noteworthy, all the above- 
mentioned estimates for bounded C 1,1 domains have the same form as for the 
ball (in this connection compare [2S1 Corollary 1.2] with [UJ Corollary 3]; see 
also © below). We also like to note that the right circular cones are merely 
special Lipschitz domains, but a number of techniques and explicit formulas 
make them an interesting and important test case (see [5J [33J 133 1231 [22] ) ■ We 
hope to encourage a further study of Lipschitz and more general domains for 
stable and other jump- type processes [2FJ Ell [3B1 [12] . We should emphasize that 
generally the estimates for Lipschitz domains cannot be as explicit as those for 
C > domains. For instance, the decay rate of harmonic and parabolic functions 
in the vertex of a cone delicately depends on the aperture of the cone, see 
[21 [51] (see also [H]). Nevertheless Lipschitz domains offer a natural setting 
for studying the boundary behavior of the Green function and the heat kernel 
of the Brownian motion and a-stable Levy processes (0 < a < 2). This is so 
because of scaling, the rich range of asymptotic behaviors depending on the 
local geometry of the domain's boundary, connections to the boundary Harnack 
principle, approximate factorization of the Green function, and applications 
in the perturbation theory of generators, in particular via the 3G Theorem, 
[Hl[l[MlE[Ii[H[Hmm[IE],and3P Theorem [TS]. Noteworthy, © is an 
approximate factorization of the heat kernel (see [TTJ [TO] in this connection). 

Cones are also examples of unbounded domains, which are only partially 
resolved by the results of [2SM27] (note that is valid only for bounded times). 
We should note that the upper bound in ([9]) was proved in 2006 by Siudeja 
for semibounded convex domains [57l Theorem 1.6] (stated for general convex 
domain in [57l Remark 1.7]). It appears that the impulse for the proof of ([9]) was 
given by Siudeja and Kulczycki in [48] Theorem 4.2], see also [H Proposition 
2.9] by Kulczycki and Banuelos. A similar but weaker upper bound was earlier 
given in (26)], see also [EH US [5j] . We also remark that 0D| Theorem 4.4] 
gives a sharp explicit estimate for the survival probability of the relativistic 
process in a half-line. Generally, the subject is far from exhausted-and it seems 
manageable with the existing techniques. 

For completeness we like to mention recent estimates [2"51 [2"T1 1531 [BUI I5"B1 [371 IB] 
for transition density and potential kernel of jump-type processes. We need 
to point out that generally these are estimates for processes without killing. 
Killing is a dramatic "perturbation" analogous to Schrodinger perturbations 
with singular negative potentials [HI [TBI EH 02] , and it strongly influences the 
asymptotics of the transition density and Green function. The asymptotics is 
crucial for solving the Dirichlet problem for the corresponding operators, see 
also [39j HO] . As we shall see, the heat kernel of the fractional Laplacian in the 
right circular cones has a power-type asymptotics at infinity, and it decays like 
the distance to the boundary to the power a/2 except at the vertex, where it 
decays with the rate of (3 € (0, a). 

The paper is composed as follows. Below in this section we recall basic 
facts about the transition density of the a-stable Levy processes killed when 
first leaving a domain. In Section [2] we give a sharp explicit estimate for the 
survival probability P x (td > t) for C > domains D. In Section [3] we prove our 
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main estimates, Theorem Q] and ([23"]). by using the ideas and results of [35J and 
[2]. Our general references to the boundary potential theory of the fractional 
Laplacian are [13] and |19j . We also refer the reader to [TS] for a broad non- 
technical overview of the goals and methods of the theory. 

In what follows, R d denotes the Euclidean space of dimension d > 1, dx is 
the Lebesgue measure on R d , and < a < 2. For t > we let pt be the smooth 
real-valued function on R d with the following Fourier transform, 

p t {x)e ix < dx = e~ m " , £,eR d . (2) 

For instance, a = 1 yields 

ft (x) = r ( (d + D/2K- (rf+1)/2 (|x|2 + ; 2)(d+1)/2 , 

the Cauchy convolution semigroup of functions [55] ■ We generally have that 

p t (x) =t- d/a Pl (t- 1/a x), x&R d ,t>0. (3) 

This follows from @. The semigroup Ptf(x) = J Md f(y)pt(y — x)dy has A Q / 2 
as infinitesimal generator ([7], [51], [13], [H]); where 

^ V ' 7r^ 2 |r(-a/2)| eio7 {|2/|>£} |y|«H-a 

Here <fi € C^ >0 (R' i ), i.e. 4> : M. d — > R is smooth and compactly supported on Mr. 
Put differently, 



p u _ s (z-i) d u 4>(u, z) + A"/ 2 i/>(u, z) dzdu — ~(f>(s,x) , 



where s S R, x S R d , and S C C °°(R x R d ) ([16 ). We denote 

the density function of the Levy measure of the semigroup {Pt} [55 ] l20 [ [15] . 
There is a constant c such that (see [20] or [9]) 

c " 1 (]^ Ard/Q )^^^ c (]^ At " d/a ) ' ^ Kd ^>°- W 

|4| and similar sharp estimates (i.e. such that the lower and upper bounds are 
comparable) will be abbreviated as follows: 

Pt (x) « A , x G R d ,t>0. (5) 
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The standard isotropic a-stable Levy process (X t , P x ) on M. d may be constructed 
by specifying the following time-homogeneous transition probability: 

P t {x,A) = / p t (y-x)dy, t>0, x £ M d , AcM d , 

J A 

and stipulating that P x (X(0) = x) = 1. Thus, P x , E x denote the distribution 
and expectation for the process starting from x. The distribution of the process 
is concentrated on right continuous functions: [0, oo) — > M. d with left limits, and 
for all s > 0, x £ M. d we have that P X (X S — X s -) = 1. It is well-known that 
(X t ,P x ) is strong Markov with respect to the so-called standard filtration [51 ITU], 
The Levy system (see [351 VII.68], [HI Appendix A], see also [S3 Theorem 2.4], 
[48l Corollary 2.8] and [3j Lemma 1]) for (X t ,P x ) amounts to the equality, 



J2f(s,Xs-,X s ) 



3<T 



= E x 



f(s,X s ,y)v(w - X s )dw ds 



(6) 



where x £ R d , / > is a Borel function on K x R d x R d , such that f(s, z, w) = 
if z = w, and T is a stopping time with respect to the filtration of X . 
For open D c K d we let 7£> = inf{f > : X t ^ D}, and we define 

V?{x,y) =pt(x,y) - E x [t d < t; p t -r D ( X T D ,y)}, x,y £ R d , t > 0, 



see, e.g., [26j[T3]. Clearly, 



< Vt{x,y) < Pt(y - a;) 



(7) 



By the strong Markov property, pf is the transition density of the isotropic 
stable process killed on leaving D, meaning that p D satisfies the Chapman- 
Kolmogorov equation: 



pf (x,z)p?(z,y)dz=p? +t (x,y), x, y £ R d , s, t > , 



and for every x £ M d , t > and bounded Borel function /, 
f(y)p?(x,y)dy = E x [T D <t; f(X t )}. 
Furthermore, for s £ R, x £ R rf , and <p £ C C °°(M x D), we have 



fi_.{x,z) \d u <f,(u t z) + ^ 2 <f>(u,z) 



dzdu 



(s,x) 



s D 



which justifies calling p D the heat kernel of the fractional Laplacian on D. In 
analogy with ([3J we have the following scaling property 



pf{x,y)=t- d ' a p t ~ 1/an {t- 1 / a x,t- 1 / a y), x,y£R d ,t>0. 



(8) 
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2 C 1 ' 1 domains 

Let D C M. d be a C 1 ' 1 domain, meaning that £) is open and there is ro > 
such that for every z € dD there exist balls B z (r ) C £> and B' z (r ) C -D c of 
radius ro, tangent at z. Denote <5d(x) = dist(x,-D c ), the distance to £> c . The 
transition density of the stable Levy process killed off D satisfies ([26 ) 



p?( X ,y)*U A 6 J>^\ U^&^^y), 0<t<l,x,yE 



(9) 



Corollary 1. If D is a C 1 ' 1 domain then 

x/2 



Px (td >t)« 1A-^U, 0<Kl,i,i/el d . (10) 



Proof. We have 



By®, 



^(tu >i) = / pf{x,y)dy. 



P x (T D >t) = I pf(x,y)dy^ lA-^-P J t (z), 0<<<l,x,ye 



where 



Clearly /t(x) < j Rd pt(x,y)dy = 1. This yields the upper bound in (|T0j) . To 
prove the lower bound we consider < t < 1 and we will first assume that 
S D (x) > i 1/Q . If \y-x\ < t 1 / a /2, then p t (x, y) w i~ d/t \ and we get 

/ t (i) > c / r d/a dy = c> 0. 

If #d(x) < t 1 /", then let z € 9.D be such that \x — z\ = Sjj(x), and consider the 
inner tangent ball B z {t x l a A ro) for D at z, with center at, say, w. We have 

f (t 1/a - \y-w\) a/2 , \ , 

It (x) > / i ^ V pt (x, y) dy . 

iB J (f 1 /«Ar ) \t 

Since pt{x, y) — i~ d//a Pi (73770 +573* )> by changing variable v = > we S e t 

It{x)> [ (l-\v\) a / 2 Pl (u,v)dv, 

JB(0,lAr ) 

where u = t~ x l a {x — w) £ B(0, 1). The latter integral is continuous and strictly 
positive for u G 5(0, 1). Thus, mi x£ £> I t (x) > 0. The proof of (fTU)l is complete. 

□ 
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Corollary 2. If D is a C 1 ' 1 domain then 

p?(x,y) ^ P x (t d > t)P y (r D > t) Pt (x,y) , < t < 1 , x, y £ R d . 

3 Cones 

For x £ R d \ {0} we denote by 6?(:r) the angle between x and the point 
(0, ...,0,1) £ R d . We fix < 9 < ir and consider the right circular cone 
r = {x £ R d \ {0} : 6{x) < 6}. Clearly, rT = T for every r > 0. By ©, 

pfc, y) = r^pfCt-^a, ^ 1/Q y) , x, y e M d , t > . (11) 

We fix xo £ r and consider the Martin kernel M for T with the pole at infinity, 
so normalized that M(x ) = 1. It is known that there is < (3 < a such that 

M{x) = \xfM(x/\x\), x^0, 

see [2J, [5T], [T7]. Since the boundary of T is smooth except at the origin, by 
[5TI Lemma 3.3], 

Af^^r^r/ 2 !^-"/ 2 , x£R d . (12) 
The following result strengthens [5J Lemma 4.2]. 
Lemma 3. If T is a right circular cone then 

P x {rr>t)^ te /2 (rV^)Al) (|i- 1/Q a;| Al)' 3_Q/2 , z e R d , t > . (13) 

Proof. Since P^ir > t) = -Pt-i/a^T > 1), we only need to prove that 

Px(tt > 1) « (<5r /2 (a;) A l) (H A l)^ 72 , x £ R d . (14) 

If |x| < 1 then HU) is a consequence of (JT2) and [2J Lemma 4.2]. If |x| > 1 
then P x {tt > 1) ~ <5r(a;) Q ^ 2 A 1. Indeed, considering C 1 ' 1 domains V and T" 
such that T' C T C T" and T" \ T' C 5(0, 1/2), we see that 5 r ,(a;) < <5 r (x) < 
(5 r »(x) < 2(5r'(x) for such x. Since P X (7T' > 1) < P x (r r > 1) < P x {tt» > 1), by 
using flT0J) we obtain (|T3|), □ 

An interesting, if trivial, consequence of (fl3|) is that 

P K (Tr > t) « P x {t t > t/2) , i > , x £ R d . (15) 

Theorem 1. 

p T t {x,y)^P x {T T >t)P y {T T >t)p t {x 1 y), x,y£R d , t>0. (16) 
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Proof. We note that the right hand side, say Rt(x, y), of (fTB]) satisfies 

R t (x,y) = t- d l a R 1 {t- 1 l a x,r 1 l a y). 

Thus, in view of (fTTj) . we only need to prove (fll))) for t = 1. 
Let F' and T" be as in the proof of Lemma [31 Then 

Pi <Pi(x,y) <Pi (x,y)- 

By © we have, for |x|, \y\ > 1, 

Pi' 0, V) « (l A (^r 72 ^)) (i A ^r /2 (y)) Pi fa, 2/) ~ Pi" 0> 2/) ■ 

Hence by Lemma [3] we obtain, 

p{(x,y) n P x (t t > l)P y (Tr > l)pi(x,y) , \x\,\y\>l. (17) 

In particular, there is a constant c such that 

p[(x,y)<cP x (Tr>l)pi(x,y), \x\,\y\>l. (18) 

If \x\ < 1 and \y\ < 4, then \x - y\ < 5 and pi(x, y) > c(l A |x - y| _d_Q; ) > c. 
By the semigroup property (J7J) and (fl"5)) . 

PiC^y) = J Pi/2( x ' w )Pi/2( w iy) dw < c J p\ /2 {x,w)dw 

= cP x {t t > 1/2) <cP x (-r r > 1) 

< cP K (r r > l)pi(x,y), |a:| < 1 , ||/| < 4. (19) 

We next assume that |x| < 1 and |y| > 4. Then |a; — y\ > 3 and pi(x,y) w 
| a; _ y \- d ~ a . Denote r x = T n 5(0, 2), r 2 = (r \ Ti) n B(0, (\y\ + l)/2) and 
T3 = r \ -B(0, (|j/| + l)/2). Using the strong Markov property and the Levy 
system © with f(s,z,w) — l^ 1 (z)lr^(w)pi_ s (w ) y) and T = 1 Arr i; we obtain 

Px(x,y) = E x [t Ti < l:,p{_ Tri (X rri ,y)} 

= p v s 1 {x,z) I v(w — z)pi_ s (w,y) dwdzds 

Jo Jr! Jr\ri 

= I I P^ixjz) / v(w — z)pi_ s (w,y) dwdzds 
Jo Jri Jr 2 

+ / P^i^tZ) v{w - z)p\_ s {w 1 y) dwdzds 
Jo JVi Jr 3 
= I + 11. 

We note that for w G I^, 

\w-y\> \y\-\w\ > |y|/4> \x-y\/8. 
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Since p^L s (w — y) < pi- s (w — y) < c(l — s)\w — y\ d we obtain 
I I I P T ^{ x i z ) I v(w — z)c; rv- — dwdzds 

II'- I *■ ' \ W - y\ d - " 



r! Jr 

-a 

-d—a I I „ri 



d+a 

< c\x — y\~ d ~ a I I p l s 1 (x,z) I v(w — z) dwdzds 
Jo Jr! Jr\r! 

= c\x-y\- d - a p x (x TTi e ryr^rr, < 1) 

< c\x-y\~ d - a P x (X TTi eT\T 1 )^M(x)p 1 {x,y). 
In the last line we used BHP ([!])• For z £ Ti and w G T^, we have 

\w - z\ > \w\ - 1*1 > \y\/2 - 3/2 > |y|/8 > \x - y|/I6, 
hence v(w — z) < c\x — y\~ d ~ a , and so 

II < c \x-y\- d - a f f ft{x,z) f p\_ s {w,y) dwdzds 
Jo Jr x Jr 3 



f f p r s 1 (x,z)dzds<c\x-y\- d - a E x Tr 1 

JO JTi 



< c\x-y\- d -< 

< cp 1 (x,y)M(x), 

where the last inequality follows from ([5]) and [2J Lemma 4.6]. By Lemma 02 

p\{x,y) < cP x {t t > l)pi{x,y), \x\ < 1, \y\ > 4. (20) 
Combining (IS]), and flU, we get 

Pi (a, J/) < cP x (t t > l)pi(x,y), x,y G K d . 
By the symmetry, semigroup property and (|15[) we obtain 

rf ( ,,,) = / /f/2 (s,»)rf /2 (», s )*» 

2 2d/Q p[(x2 1/Q , w2 1/a )p{(w2 1/a 1 y2 1/a )dw 



< cP x2 i/c{t v > l)P y2 i /Q (r r > 1) / p 1 / 2 (x,w)p 1 /2(w,y)dw 

< cP x (r r > l)P y (rr > l)pi(x,y). 

We will now prove the lower bound in (fl6|) . We first assume that |a;| < 1 
and \y\ < 2, and we let T4 = F n B(0, 4). Since T4 is bounded, the semigroup 
p^ 4 is intrinsically ultracontractive ([47]). In particular, 

Furthermore, by J5[ Lemma 4.6] and Lemma[3]we obtain 

B*tt 4 >cM(y) >cP y (r r > 1) . 
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We see that 

p\ /2 {x,y) > cP x (r r > 1)P v {tt > l)pi{x,y), \x\ < 1, \y\ < 2. (21) 
If \x\ < 1 and \y\ > 2, then by the semigroup property, (fT7|) and ([2"Tj) . 

r ^ ^jr 

r 

2„ 



p 1 (x,y) = j p 1/2 {x,z)p 1/2 (z 1 y)dz 

> cP x (t t > l)P y (T T > 1) / P z (t t > l) 2 pi{x,z) Pl {z,y)dz 

Jri\B(p,i) 



> cP x {t t > l)P y (Tr > l)pi(x,y) [ P z (T r >l) 2 dz. 

^ri\s(o,i) 



Hence 



p r 1 (x,y)>cP x (T r >l)P y (T r >l)p 1 (x,y), \x\ < 1, \y\ > 2. (22) 
By (TT7)) . (|2T|) . (|2^|) , symmetry (and scaling), we get the lower bound in (TIB)) . □ 
We like to note that Theorem [T] strengthens [2 Corollary 4.8]. Also, 

p T t {x,y)Kp T t/2 (x,y), x,y € R d , t > 0. 
In view of Lemma [3J for the right circular cone T, (|16p is equivalent to 

rfCx.^^ft-^Al) (l^^l Al)^ /2 (t- d/a A |a; _* |d+a )(23) 

Pr 1 /"!/) A l) {\t- 1 ' a y\ A l)^" 72 , i>0,x je i d . 



This is explicit except for the exponent (3 (see [5] in this connection). Recall 
that J °° p^(x,y)dt = Gr(x,y), the Green function of T. By integrating (|23j) 
one can obtain sharp estimates for the Green function of the right circular cone. 
For d > 2 the estimates-first given in [5TJ Theorem 3.10]-are the following, 

^^ 1A j<!^^y-n „, e *. (24) 

We skip the details of the integration (similar calculations are given in [57] and 
[2l>]). Noteworthy, j3 = a/2 if T is a half-space [2J. For the case of dimension 
d= 1, and T = (0,co), we refer the reader to [24], see also [26] Corollary 1.2]. 

As stated in Introduction, we expect to be true quite generally. In 
particular the approximation should hold for domains above the graph of a 
Lipschitz function for all times t > 0. Corollary [2] confirms this conjecture for 
C 1 ' domains and small times, while Theorem [1] proves it for the right circular 
cones and all times. By inspecting the relevant proofs in [51] . the reader may 
also verify without difficulty that Theorem [T] and (|24[) hold the same for all 
those generalized cones (PI) in W\ d > 2, which are C 1,1 except at the origin. 
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On the other hand, if D is a bounded C 1 ' 1 domain, and if we denote by — Ai 
the first eigenvalue of A"/ 2 on D (i.e. when acting on functions vanishing off 
D), then by the intrinsic ultracontractivity (see, e.g., Theorem 1.1]), 

p?(x,y)^S a D /2 (x)S^ 2 (y)e^ t , t>l,x,yeR d , 

and so 

P x (t d >t)K 8 a J 2 {x)e- Xlt , t>l,xeR d . 

Therefore |T]) fails for large times t if D is bounded (see also [47]). 
Acknowledgments. Wc thank Michal Ryznar and Mateusz Kwasnicki for 
discussions and suggestions. 
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